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Introduction 

The fundamental role in the theory of linear elliptic equations and systems is played 
by results on regularity of solutions up to the boundary of a domain. The following 
classical example serves as an illustration. 
Consider the Dirichlet problem 

[ Au = in O, 

\ Trit = / on dn, 

where f2 is a bounded domain with smooth boundary and Tr u stands for the boundary 
value (trace) of u. Let u S Lp(Q), 1 < p < oo, that is 

\Vu\ p dx < oo 

n 

and let / belong to the Besov space Bp~ 1 ^ p (dfl) with the seminorm 



dQ Jdtl 



\x - y|«+P- 2 x v 
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It is well-known that 

TrL*(Sl)=B$-V p (d[i). 

Moreover, the harmonic extension u of / G Bp 1 ^ p {dVL) belongs to L p (fl) and the 
norm of the gradient Vit in L p (f2) is equivalent to the above seminorm in Bl- 1/P (dn). 
This fact highlights the following topics of interest: 

• Replace the domain with smooth boundary by a domain in more general class 
and study the effect of irregularities of dft. 

• Replace A by a more general elliptic operator with variable coefficients, and 
study the impact of low regularity assumptions on the coefficients. 

• Understand the correlation between the smoothness of data and the smoothness 
of solutions. 

Another theme, somehow related to these, is: 

• Describe the local behaviour of solutions near a boundary or interior point of 
the domain fl. 

In the present paper we survey results in the directions just mentioned. Most of 
them were obtained during the last decade and concern elliptic equations and systems 
of arbitrary order. Special attention is paid to the Stokes system. The selection of 
topics is partly influenced by our involvement in their study. Here is the plan of the 
article. 

Section 1 is dedicated to weak solutions with Besov boundary data, with coeffi- 
cients of differential operators and the unit normal to the boundary in classes close to 
BMO. Results of a similar nature for the Stokes system are discussed in Section 2. 
In Section 3, strong solutions in Sobolev spaces are considered. Here sharp additional 
conditions on the Lipschitz boundary are reviewed. In particular, in Subsection 3.2 
we speak about strong solvability of the Stokes system. In Section 4, asymptotic 
formulas for solutions of the Dirichlet problem near an isolated point of the Lipschitz 
boundary and at a point in the domain are dealt with. 

1 Weak solvability of higher order elliptic systems 
with coefficients close to BMO in Lipschitz do- 
mains 

1.1 B ackground 

The present section is mostly based on results of the paper by V. Maz'ya, M. Mitrea 
and T. Shaposhnikova M MSlj . We start with mentioning earlier works. 

The basic case of the Laplacian in arbitrary Lipschitz graph domains in R", is 
treated in the work of B.Dahlberg and C.Kenig |DK| . in the case of L p -data, and 
D. Jerison and C. Kenig |JKj . in the case of data in the Besov space B^ with < s < I. 
The local regularity in the Sobolev class Wp of solutions to second order equations 
with coefficients in VMODLoo was established by F. Chiarenza, M. Frasca, P. Longo 
[CFL] , 
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In spite of substantial progress in recent years, there remain many basic open 
questions for higher order equations, even in the case of constant coefficient operators 
in Lipschitz domains. One significant problem is to determine the sharp range of p's 
for which the Dirichlet problem for strongly elliptic systems with L p -boundary data 
is well-posed. In |P V) . J. Pipher and G. Verchota have developed a L p -theory for real, 
constant coefficient, higher order systems 

\a\—2m 

when p is near 2, i.e., 2 — e < p < 2 + e with s > depending on the Lipschitz 
character of fl. On p. 2 of [PVj the authors ask whether the L p -Dirichlet problem for 

these operators is solvable in a given Lipschitz domain for pe (2 — e, 2 ^Zp + £ ): an d 
a positive answer has been given by Z. Shen in |Shj . Let us also mention the work 
|AP] of V. Adolfsson and J. Pipher who have dealt with the Dirichlet problem for the 
biharmonic operator in arbitrary Lipschitz domains and with data in Besov spaces, as 
well as [Vej where G. Verchota formulates and solves a Neumann-type problem for the 
bi-Laplacian in Lipschitz domains and with boundary data in L2. In [MMT D. Mitrea, 
M. Mitrea, and M. Taylor, treat the Dirichlet problem for strongly elliptic systems 
of second order in an arbitrary Lipschitz subdomain of a (smooth) Ricmannian 
manifold and with boundary data in Bp(dQ), when 2 — e < p < 2 + e and < s < I. 

We mention some recent results on the second order elliptic equations and systems 
with coefficients in VMO due to G. Di Fazio [DP], L. Caffarelli and I. Peral [CP], B. 
Stroffolini [St], D. Guidetti |Guj . for higher order equations see also D. Palagachev 
and L. Softova [PS) . 

1.2 Domains and function spaces 

Let us turn to the article [MMS1] . We make no notational distinction between spaces 
of scalar-valued functions and their natural counterparts for vector-valued functions. 

Recall that a domain f2 is called Lipschitz graph if its boundary can be locally 
described by means of (appropriately rotated and translated) graphs of real- valued 
Lipschitz functions. 

It is shown by S. Hofmann, M. Mitrea, and M. Taylor in jHMTj that fi is a 
Lipschitz graph domain if and only if it has finite perimeter in the sense of De Giorgi 
(see |DGlj . [Fe] . |BMj ) and (i) there are continuous (or, equivalently, smooth) vector 
fields that are transversal to the boundary and (ii) the necessary condition dfl — dVl 
is fulfilled. 

By |HMT| , a bounded nonempty domain of finite perimeter for which dfl = dfl is 
Lipschitz graph if and only if 

inf{||i>-w|k (8n) : w = (wi, . . . , w„) e C°(dn), \u\ = 1 on dfl} < 

with v being the outward normal to dfl. 

Everywhere in this section we assume that f2 is a Lipschitz graph domain in R™, 
with compact closure f2 and with outward unit normal v — (v\, ...,v n ). Let m be an 
integer. Consider the operator 

£(X,Dx)U:= D a {A a p{X)D f3 U), (1.1) 

|o|=|(8|=m 
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with the data ^ 

^i=9k on ^ 0<*<m-l. (1.2) 

The coefficients A a p are square matrices with measurable, complex-valued entries, for 
which 3 k > 

£ IIA^IU^n) <« _1 (1.3) 

|a| — 1/3|— m 

and such that the coercivity condition 



Re / V (A a0 (X)D^V(X),D a V(X))dX>K V 



2 

Il 2 (o) 

/f2 | a | = |/3|= m | Q | 



holds for all complex vector- valued functions V G C^°(f2). 

Let W belong to the usual Sobolev space WZ™(fi). It is natural to take 

V^ a Tr[TO], 0<fc<m-l, 

\a\—k 

where v a := is" 1 • ■ • if a — (ai, a n ). Now, if p € (1, oo), a e (— 1/p, 1 — 1/p) 
are fixed and 

p{X) := dist (X,dft) 
the space W^ n ' a (n) is defined as the space of vector- valued functions for which 

V f \D a U(X)\ p p(Xy a dX < oo. (1.4) 

Further we set 

V™< a (n) := the closure of C °°(O) in W;"< a (0) (1.5) 
and introduce the dual space 

V- m "(fl) := (^; i - a (^))* 

For any U € W™' a (0), the traces of D a l4, < |a| < m- 1, exist in B*(dSl), where 
s := 1 — a — 1/p, < s < 1 (see E. Gagliardo |Gag| for a — and S. Uspenskh |Usp| ). 
Recall that / e L p (dil) belongs to Bp(dft) if and only if 

r i/w-/(y)r ^ 

The above definition takes advantage of the Lipschitz manifold structure of dfl 
which allows one to define smoothness spaces of index less than 1. The approach is 
no longer effective when the order of smoothness exceeds 1. 

Let us illustrate the necessity of working with boundary data different from those 
in spaces of traces of usual Sobolev spaces by considering Dirichlet problem for the 
biharmonic operator 

WGW|(0), A 2 U = 0mn, 

TvU = g on 90, (v, Tr [VW]) = gi on dfl. 
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One might be tempted to believe that a natural class of boundary data is B 2 (dfl) x 
B^ 2 (dfl), where by definition B^ 2 {dfl) and B^ 2 (dfl) are the spaces of traces of 
functions in W|(f2) and W^fii), respectively. 

However, this formulation has serious drawbacks. The first one is that the mapping 
W|(0) 3W^(!/, Tr [VW]) € B 2 1/2 (af7) 

1 /2 

is generally unbounded. In fact, its continuity implies v £ B 2 (dQ) which is not 
necessarily the case for a Lipschitz domain, even for the square S = [0, l] 2 . 

Second, this problem may fail to have a solution when (<7o> <?i) is an arbitrary pair 
in B^ 2 (d£l) x B^ 2 (dfl) . Indeed, consider the case il = S and g = 0, g\ = 1. It is 
standard that the main term of the asymptotics near the origin of any solution U in 
W 2 {S) is given in polar coordinates (r, 9) by 

-^((^f)sin^costf). 

Since this function does not belong to W 2 (S), there is no solution in this space. 

A new point of view has been introduced by H. Whitney in |Whj who considered 
higher order Lipschitz spaces on arbitrary closed sets. An extension of this circle of 
ideas pertaining to the full scale of Besov and Sobolev spaces on irregular subsets of R" 
can be found in the book |JWj by A. Jonsson and H. Wallin. The authors of [MMSlj 
further refined this theory in the context of Lipschitz domains. The description of 
higher order Besov spaces on the boundary of a Lipschitz domain il C M™ in [MMS1 
runs as follows. 

For m £ N, p £ (l,oo), s £ (0,1), the space B™~ 1+8 (dQ) is introduced as the 
closure of 

{(D a V\ dn )\ a \< m -i ■■ VeC °°(R")} 

in Bp(dil). An equivalent characterization of B™~ 1+s (dil) which involves higher 
order Taylor remainder in place of f(X) — f(Y) in (|1.6p runs as follows (see Sect. 7.1 
of [MMSlp . 

For a collection of families / = {f a }\ a \<m—i of measurable functions defined on 
dfl, there is the equivalence relation 

II/Hb™-^ 8 ^) ~ ^2 Wf<x\\L p (dCl) (1-7) 

|a| <m— 1 

/ f f \R a {X,Y)\P , , \Vp 

+ > / / 1 1 7 f; u i Tda x da Y ) < oo, 

, ^ ^JanJan I - y pK 1 +-H)+-i x y ) 

\oc\<m— 1 

where 

R a (X,Y) :=f a (X)- ]T ^U+p(Y)(X-Yf, X,Y £ dQ, (1.8) 

\/3\<m-l-\a\ P ' 

It is standard to prove that £?™~ 1+s ((9f2) is a Banach space. Also, trivially, for any 
constant n > 0, 

H/a|U,(flO)+ (// \ X - yS*i?i'+«-|a|)+n-l ^X^y) ' (1.9) 

|a|<ro-l |a|<m— 1 x.veon 

\X-Y\<k 
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is an equivalent norm on B™ 1+s (dft). 

In order to formulate the trace and extension theorem for the spaces B™~ 1+s (dft) , 
we first give its analogue for lower smoothness spaces which is essentially due to 
S.Uspenskii |Usp| . 

Lemma 1 For each 1 < p < oo, —1/p < a < 1 — 1/p and s := 1 — a — 1/p, the trace 
operator 

Tr : W^ a (Q) — > B s p (dQ) (1.10) 

is well-defined, linear, bounded, onto and has Vp< a (£l) as its null-space. Furthermore, 
there exists a linear, continuous mapping 

£ : — ► W*' a (fl), (1.11) 

called extension operator, such that Tro£ = / (i.e., a bounded, linear right-inverse of 
trace). 



For higher smoothness see the following assertion which is Proposition 7.3 in 
jMMSlj . 

Proposition 1 For 1 < p < oo, —l/p<a< 1 — 1/p, s := 1 — a — 1/p € (0, 1) and 

m G N, define the higher order trace operator 

tr m _! : W™> a (Q) — > B™- 1+S (dn) (1.12) 

6t/ setting 

tr m ^U:= \i^Ti[D a U]\ , (1.13) 

l J |a|<m — 1 

where the traces in the right-hand side are taken in the sense of Lemma [7J Then 
(| 1 . 1 2|) is a well-defined, linear, bounded operator, which is onto and has V™' a (£l) as 
its null-space. Moreover, it has a bounded, linear right-inverse, i.e., there exists a 
linear, continuous operator 

e ■ B™- 1+S (dn) — > w™> a (ri) (i.i4) 

such that 

f = {fa}\ a \< m -i G B™- 1+S (dn) »M Tr [!>"(£ /)] - /„. (1.15) 

Now, a necessary condition for the boundary data {<?fc}o<fc<m-i in the Dirichlct 
problem (|1.2|) is that 

3 / = {/a}|a|<m-l G S^+^df!) Such that 

fc! (1.16) 
gk = y — r 1/0 /a for each < k < m — 1. 

Q | = A: 

This family {<7fc} is organized as a Banach space, and is denoted by W™~ 1+s (dH). 

The space takes a particularly simple form when m — 2. To describe it we need 
the notation for the tangential derivative d/drjk given by 

d d d 

»3ltZ7 ~ "kjr-, l<j,k<n (1.17) 



drjk ' dx k dxj 



G 



and the tangential gradient on the surface dQ 

Vtan := Vjd/dTjk)l<k<n- 
3 

Then, for each Lipschitz graph domain f2 C 1" and each 1 < p < oo, s G (0, 1), 

W^ +S (dn) = {(go,gi) e Ll(dSl) ® L p (dty : vgi+Vt^ga € ^(SQ)}. (1.18) 

This has been conjectured to hold (when s = 1 — 1/p) by A. Buffa and G. Geymonat 
on p. 703 of [BG]. 

1.3 Formulation of the Dirichlet problem 

Broadly speaking, there are two types of questions pertaining to the well-posedness of 
the Dirichlet problem in a Lipschitz domain Q for a divergence form strongly elliptic 
system of order 2m with boundary data in W£ l ~ 1+s (dQ,). 

Question I. Granted that the coefficients of £ exhibit a certain amount of smoothness, 
identify the indices p, s for which this boundary value problem is well-posed. 

Question II. Alternatively, having fixed the indices s and p, characterize the smooth- 
ness of 9f2 and of the coefficients of £ for which the aforementioned problem is well- 
posed. 

Both questions are discussed in |MMSlj . 

B. Dahlberg (see [Dahlj has constructed a domain fl c K™, bounded with C 1 - 
boundary, and a function / e C 00 (fj) such that, for each p £ (1, oo), 

Au = f,uG W£(tt), u\ dn = =>■ d 3 d k u i L p (n), 

where dj = d/dxj. In other words, Dahlberg's domain ft has the property that for 
each p G (l,oo), the Poisson problem for the Laplacian with homogeneous Dirichlet 
boundary conditions fails to be well posed, in the sense that it is no longer reasonable 
to expect solutions with two derivatives in L p (f2). 

A fundamental issue raised by this counterexample: 

Identify those Sobolev-Besov spaces within which the natural correlation between 
the smoothness of the data and that of the solutions of the Dirichlet problem is pre- 
served when the domain in question is allowed to have a minimally smooth boundary. 

Consider the Dirichlet problem for the operator (|1.1[) 

r C(X,D X )U — T in ft, 
i d k lA 

— — r- = Ofc on dQ, < k < m — 1. 

Proposition 2 J/We W£ l ' a (Q) then, necessarily, 

T e V- m ' a (Q), g := {g k } < k < m -i G W™- 1+S (dn) 

and, moreover, 

\\9\\w?- 1 +'(da) + W^Wv^iu) ^ c \\ u \\w™- a (n)- 
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For the results in the converse direction, the main hypothesis in |MMSlj requires 
that, at small scales, the so called local mean oscillations of the unit normal to dCl 
and of the coefficients of the differential operator C{X, Dx) are not too large, relative 
to the Lipschitz constant of the domain CI, the ellipticity constant of C{X,Dx), and 
the indices of the corresponding Besov space. 



By the local mean oscillation of F E we understand 



{F} 



Osc(O) 



:= lim ( 



JB^nn J B r 



sup 

{B e } n JB E nn JB e nft 



F(x) - F(y) 



dxdy ) , 



where {-B e }fi stands for the family of balls of radius e centered at points of ft and 
the barred integrals denote the mean values. Similarly, the local mean oscillation of 
/ G Li(0fi) is 



{/} 



osc(an) 



lim ( sup + 



B E ndn JB E ndn 



m - f(y) 



ds x ds y ^j , 



where {B e }dn is the collection of n-dimensional balls of radius e with centers on dCl. 

Note that smallness of the local mean oscillation {v}o sc (dn) does not imply small- 
ness of the Lipschitz constant. Indeed, let 



where 



Then ||^|| Lo 



n = {(x,y)eR 2 ,y>^ £ (x)}, 

ip e {x) — x sin(e log | | 1 ) . 
1, while \\<p' £ \\bmo(r) < C e. 



1.4 Solvability of the Dirichlet problem in W™ ,a 

The main result in M.MSTj runs as follows. 

Let fl C M. n be a bounded Lipschitz domain whose Lipschitz constant is does not 
exceed M, and assume that the operator C{X,Dx) of order 2m is strongly elliptic, 
and has bounded, measurable complex coefficients. 

Theorem 1 There exists a positive constant C , depending only on M and the ellip- 
ticity constant of C, such that: For each p £ (1, oo) , s £ (0, 1) and a := 1 — s — 1/p, 
the Dirichlet problem 

( D a {A a p{X)D f *U)=J : for X eil, 

1 (1.19) 
d k Li 

——77 = gk on dft, < k < m — 1. 
ov K 

with T € V~ m ' a (n) and g := {gk}o<k<m-i in W™~ 1+ '(d£l) has a unique solution 
U G W™' a (Cl) if the coefficient matrices A a p and the exterior normal vector v to dfl 
satisfy 

Mo + { A "/3>o S c(n) (1-20) 

\a\ — \/3\—m 



<Cs(l- s)(p 2 {p - l)- 1 + «-*(! - s)" 1 ) \ 
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For second order operators the factor s(l — s) in the last inequality can be removed. 
Furthermore, there exists C — C(dVL, A,p, s) > such that 

\M\w™'* m < C(\\g\\ w ™-^ {dn) + m v -^ m )- (1-21) 
The next assertion, obtained in jMMSlj , is a byproduct of the proof of Theorem 

m 

Theorem 2 Let O C R n be a bounded Lipschitz domain whose Lipschitz constant 
does not exeed M, and assume that the operator C(X,Dx) of order 2m is strongly 
elliptic, and has bounded, measurable (complex) coefficients. Then there exists e, 
depending only on M and the ellipticity constant of C, such that: 

For each p G (1, oo), s G (0, 1) and a := 1 — s — l/p, the Dirichlet problem for C 
with J- G V~ m >°(fi) and g := {.9£:}o<fc<m-i in W™ +s {d£l) has a unique solution 
U G \\','," "M) if 

I2" 1 - p~ x \ < e and \a\ < e, 

Recently M. Agranovich |Agl| obtained this type of results (\p — 2| is small and 
< s < 1) for both Dirichlet and Neumann problems for a subclass of strongly elliptic 
systems with Douglis-Nirenberg structure in bounded Lipschitz domains. The regu- 
larity results in |Agl| concern solutions in spaces of Bessel potentials and Besov 
spaces Bp with coefficients of differential operators satisfying the uniform Lipschitz 
condition. The approach is based on regularity methods due to Savare as well as on 
author's developement of interpolation theory of spaces Hp and Bp with a of arbi- 
trary sign (see |Ag2| for more details), where essential role is played by an extension 
operator from O to R n introduced by V. Rychkov |Ry| . To be specific, Agranovich 
considered a Douglis-Nirenberg system with the principal part Lq whose entries are 
given by 

L j>k (x,D)= \J /rf^ ;,i/r) 

The coefficients g(x) are complex- valued and the formal self-adjointness of the 
operatot Lj^ is not assumed. The principal symbol of the system, i. e. the matrix 
Lo(x,£) with entries Lj^(x,^) is subordinate to the condition of strong ellipticity: 

ReZoM > CA(0, 

where A(^) is the diagonal matrix with entries |£| 2m J on the main diagonal, j = 
1, . . . , I, and C is a positive constant. The use the Savare method requires the addi- 
tional condition 

/3|=™ fc 

for any numbers Q at all points x E fl. 

Let BMO and VMO stand, respectively, for the space of functions of bounded mean 
oscillations and the space of functions of vanishing mean oscillations (considered either 
on il, or on dfl). It can be proved that 

{F} 0sc ~dist (F.VMO) 
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where the distance is taken in BMO. Thus the small oscillation condition introduced 
in Theorem [T] holds if 

dist (u, VMO) + dist ( A <*P ' VM0 ) 

\a\ = \/3\= m 

< Cs(l - s)(p 2 (p- l)- 1 + s-^l - s)- 1 )" 1 . 

This is the case if, e.g., [v]bmo + SI^WIbmo is sufficiently small hence, trivially, if 
v e VMO(dQ) and A a p belong to VMO(fi), irrespective of p, s, C and £1 

Other examples of domains satisfying the hypotheses of Theorem Q] are: Lipschitz 
domains with a sufficiently small Lipschitz constant, relatively to the exponents p and 
s. In particular, Lipschitz polyhedral domains with dihedral angles sufficiently close, 
depending on p and s, to ir. 

The innovation in [MMS2 that allows to consider boundary data in higher-order 
Besov spaces, is the systematic use of weighted Sobolev spaces. In relation to the 
standard Besov scale, we would like to point out that 

a = 1 - s - - G (0, 1 - 1/p) W" l - a (Q) S™- 1+ " +1 /P(n), 
p 

and 

a = i-s-~e f-i/p,o) => b," 1 ~ 1+s+1/p (^) ^ w™< a (n). 

p f V 

Of course, W™'°(f2) is just the classical Sobolev space W™(f2) when a = 0. 
In the context of Theorem [TJ for a Lipschitz dfl, 

J2 \\Tr[D a U}\\ B , (dn) -(j2 fj{Xy {1 - s) - 1 \D a U{X)\vdx) 1 '\ 

|o:|<m— 1 |a|<T?i 

uniformly in U satisfying C(X,Dx)U = in fL 

This generalizes the trace and extension result by S. Uspenskh |Usp| according to 
which 




if 1 < p < oo and < s < 1 , uniformly for u harmonic in the upper-half space. 

Of course, condition (( 1 . 3|) ensures that the left-hand side of (jl.20[) is always finite 
but it is its actual size which determines whether for a given pair of indices s, p, the 
problem (|1.2[) . (|1.4j) . (|1.16|) is well-posed. Note that the maximum value that the 
right-hand side of (|1.20[) takes for < s < 1 and 1 < p < oo occurs precisely when 
p = 2 and a := 1-s-l/p = 0. As (fL~2"0"l) shows, the set of pairs (s, 1/p) € (0, l)x (0, 1) 
for which p.2j) is well-posed in the context of Theorem [1] exhausts the entire square 
(0, 1) x (0, 1) as the distance from v and the to VMO tends to zero (while the 

Lipschitz constant of S7 and the ellipticity constant of C stay bounded). That the 
geometry of the Lipschitz domain £1 intervenes in this process through a condition 
such as (|1.20p confirms a conjecture made by P. Auscher and M. Qafsaoui in |AQ| . 

Theorem [2] can be viewed as a far reaching extension of a well-known theorem of 
N.Meyers, who has treated the case m = 1, I = 1 in |Mey| . The example given in 
Section 5 of |Mey| shows that the inclusion of p into a small neighborhood of 2 is 
a necessary condition, even when dfl is smooth, if the coefficients A a p are merely 
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bounded. For higher order operators we make use of an example due to V. Maz'ya 
[Mazl (cf. also the contemporary article by E. De Giorgi |DG2] V Specifically, when 
m G N is even, consider the divergence-form equation 

A5 m - 1 £ 4 A5 m - 1 M = in fi := {X e K" : \X\ < 1}, (1.22) 

where £4 is the fourth order operator 

U{X.D X )U := aA 2 U + bjr a(&^«)+6^ 9^Saw) 

+c ± d kdl ( Xi f x f Xl d idj u). (1.23) 

i,j,k,l=l ' ' 

Obviously, the coefficients of La(X, Dx) are bounded, and if the parameters a, b, c 6 
M, a > 0, are chosen such that b 2 < ac then C along with A2 m_1 £ 4 A2 m_1 are 
strongly elliptic. Now, if W£ denotes the usual L p -based Sobolev space of order s, 
it has been observed in jMazlj that the function U(X) := \X\ e+m - 2 e W 2 m (f2) has 
TrU G C co (dft) and is a weak solution of (| 1 . 22[) for the choice 



n n 2 (n — l)(bn - 



(1.24) 



2 V4 a + 26 + c 

Thus, if 

a:=(n-2) 2 + e, b:=n(n-2), c := n 2 , e > 0, 
the strong ellipticity condition is satisfied and — 9(e) becomes 

2- ? V2 + ne 1/2 /2(4(n-l) 2 +e) 1/2 . 

However, U € W™(Q) if and only if p < n/{2 - 6(e)), and the bound n/(2 - 9(e)) 
approaches 2 when e — > 0. An analogous example can be produced when m > 1 
is odd, starting with a sixth order operator Cq(X, Dx) from [Mazlj . In the above 
context, given that 1^(0) VMO(57), it is significant to point out that both for 
the example in |Mey| , when n = 2, and for (| 1 . 22|) when n > 3, the coefficients have 
their gradients in weak-L n , yet they fail to belong to W^(fl). 

We conclude this subsection with a remark pertaining to the presence of lower order 
terms. More specifically, granted Theorem [lj a standard perturbation argument (cf., 
e.g-, [H]) proves the following. Assume that 

A(X,D X )U:= J2 D a (A a ( j (X)D f3 U), leSl, (1.25) 

O<\ot\,\0\<m 

where the principal part of A(X, Dx) satisfies the hypotheses made in Theorem[T]and 
the lower order terms are bounded. Then, assuming that (11.20)) holds, the Dirichlet 
problem (|1.19[) is Fredholm with index zero, in the sense that the operator 

w™> a (n)3U^ (a(x,d x )u , {d k u/dv k } a < k < m ^ e v~ m < a (n)@w; n ~ 1+s m) 

is so. Furthermore, the estimate 

\\U\\w^ m < C (\\n Vp -™>° { n) + \\9\\w--^(on) + Wh pm ) (1-26) 
holds for any solution U € % m - a (ft) of (fl~T9ll . 
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1.5 Comments on the proof of Theorem \T\ 

One difficulty linked with the case m > 1 arises from the way the norm in 

(J] / \D a U(X)\Pp(Xrdx) 1/P <oo. (1.27) 

|a|<m ^ 

behaves under a change of variables 

x : n = : X n > <p{X')} -> R™ 

destined to flatten the Lipschitz surface (517. When m = 1, a simple bi-Lipschitz 
changes of variables such as fl 9 (X',X„) M> (X',X n — (f(X')) G R™ will do, but 
matters are considerable more subtle in the case m > 1. The extension operator used 
in [MMSlj was introduced by J.Necas (in a different context; cf. p. 188 in [Nee] ) and 
then rediscovered by V. Maz'ya and T. Shaposhnikova in MSl] (see also [MS2;. and 
later by B. Dahlberg, C. Kenig J. Pipher, E. Stein and G. Verchota (cf. |Dah2] and the 
discussion in [DKPV ), and S. Hofmann and J. Lewis in |HL| . 

The extension operator in question is introduced in the following way. Fix a 
smooth, radial, decreasing, even, non- negative function ( in R™ -1 such that £(t) = 
for \t\ > 1 and 

C(t)dt = l. (1.28) 

(For example, ((t) :— cexp(— 1/(1 — |i| 2 )+) for a suitable c.) Define the extension 
operator T by 

(T<p)(x',x n ):= J C(t)<p(x'+x n t)dt, (x',x„)eK+, (1-29) 

acting on functions <p from Li i ; oc (R' i_1 ). 

The two estimates below provide useful properties of the operator T. 

(i) For each mufti- indices a with |a| > 1 there exists c > such that 

D°{T<p){x)\ < «!r H [V^MO(r->), Vx = (x',x n ) G R™. 

(ii) If V x ,(p £ BMO(R"- 1 ) then V(T^) G BMO(R™ ) and 

[V(T^)] BMO (k«_) < c [V :C '(p]bmo(r™- 1 )- 

Another ingredient in the proof of Theorem[T]is deriving estimates for D"D@G(x, y) 
where G is the Green matrix of the operator 

L(D X )= ]T A a D a x , 

\ol\— 2m 

i.e. a unique solution of the boundary- value problem 
L(D x )G(x, y) = S(x - y)h for x G R" , 
(— J G){{x',Q),y) = QI l forx'GR"- 1 , 0<j<m-l, 
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where y € R" is regarded as a parameter. The methods employed in earlier works 
are based on explicit representation formulas for G(x, y) and cannot be adapted to 
the case of non-symmetric, complex coefficient, higher order systems. The approach 
in [MMS1] consists of proving that the residual part R(x, y) := G(x, y) — $>(x — y), 
where $ is a fundamental matrix for L(D X ), has the property 

\\D2DPR(x,v)\\<C\x-v\-» 

for \a\ = \j3\ = m, x,y € WL, where y := (y' , —y n ) is the reflection of the point y £ 
with respect to <9R" . 

2 Stokes system 

By Bp q (M. n ) we denote the space of functions in R" having the finite norm 

Hk,(«-) = (/ IIAfcVMttll^^lfcr-^dfcJ^ + lltill^.,^) , (2.1) 

where {s} > 0, p, q > 1 and AhV = v(- + h) — v(-). For p = q we use the notation 
Bp(R n ). 

The Besov scale Bp* (f2) is defined by restricting the (tempered) distributions from 
the corresponding spaces in W 1 to the open set f2. Also, B* (dtl) stands for the Besov 
class on the Lipschitz manifold dQ, obtained by transporting (via a partition of unity 
and pull-back) the standard scale B^ g (]R n_1 ). 

2.1 Weak solvability in Besov and Triebel-Lizorkin spaces 

Consider the Stokes system in an arbitrary bounded Lipschitz domain Q C M", n > 2, 



Ait - Vtt = / e Bp + q lp 2 {n), &xvu = g<EB ! £q p 1 (CL), 

neBZ P (SI), Tru = heB s M (dty, 
subject to the (necessary) compatibility condition 

/ (p, h) da — / g(X) dX, for every component O of Q. (2-3) 
JdO Jo 



When Oft is sufficiently smooth (at least of the class C 2 ), the problem (|2.2p was 
studied in many papers, first in Sobolev spaces with an integer amount of smooth- 
ness by V.A. Solonnikov |Solj . L. Cattabriga [Cat] . R. Temam |Tem] . Y. Giga |Gig| , 
R. Dautray and J.-L. Lions [DL , among others. This has been subsequently extended 
by C. Amrouche and V. Girault |AG] to the case when 917 £ C 1,1 and, further, by 
G.P. Galdi, C.G. Simader and H. Sohr |GSS) when dfl is Lipschitz, with a small Lips- 
chitz constant. The case when Q is a polygonal domain in K 2 , or a polyhedral domain 
in R 3 also has a reach history. An extended account of this field of research can be 
found in V. Kozlov, V. Maz'ya and J. Rossmann's monograph [KMRj . which also con- 
tains references to earlier work. Among recent publications we mention the paper 
by V. Maz'ya and J. Rossmann MRl] as well as their book |MR2] . Lipschitz and C 1 
subdomains of Riemannian manifolds were treated in the paper [DMij by M. Dindos 
and M. Mitrea, and the paper [MT] by M. Mitrea and M. Taylor. 
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The principal result in MMS2 on solutions of the problem (|2.2p holds under the 
mild condition on the normal v to the boundary of the Lipschitz graph domain and 
runs as follows. 

Theorem 3 Assume that 

71-1 



<p<oo, 0<g<oo, (n- 1)(| - l) + < s < 1. (2.4) 

Then there exists S > which depends only on the Lipschitz character of and 
the exponent p, with the property that if {^}osc(ar2) < S, then the problem \2. l 2$ is 
well-posed (with uniqueness modulo locally constant functions in D, for the pressure). 
There exists a finite, positive constant C = C(Q,p,q,s,n) such that 

\\u\\ s+ i + inf ||tt — c|| s+ i_ a 
B P , q p (n) c Bp , q p (n) 

<c(\\f\\ .+i_ 2 +\\g\\ .+!_! +\\h\\ B - on)), 
V B p , q p (o) s p ,, p (n) v 

where the infimum is taken over all locally constant functions c in Q. 

Moreover, analogous well-posedness result holds on the Triebel-Lizorkin scale (for 
its definition see, for example, the book by T. Runst and W. Sickel [R5]), i.e. for the 
problem 

Am - Vtt = / e Fpp~ 2 (Q,), divu = g£ Fpp' 1 (ft), 

ueFpp(n), TieFpp (O), Tiu = heB^ q (dn), 

This time, in addition to the previous conditions imposed on the indices p, q, it is 
also assumed that p, q < oo. 

It should be noted that conditions (|2.4|) describe the largest range of indices p, q, s 
for which the Besov spaces B s v q (dCl) can be meaningfully defined on the Lipschitz 
manifold 90. 

Regarding the nature of the main result, note that no topological restrictions are 
imposed on the Lipschitz domain (in particular, the boundary can be disconnected). 
This is significant since the approach is via boundary layer potentials, whose invert- 
ibility properties are directly affected by topological nature of the domain. That 
difficulty is overcome by using certain mapping properties of the hydrostatic layer 
potentials: single layer potential S and double layer potential T>\ for the velocity and 
single layer potential Q and double layer potential V\ for the presure. We list their 
properties in the following theorem. 

Theorem 4 Let Q be a bounded Lipschitz domain in R™, n > 2, and assume that 
A e M, ^ < p < oo, (n - l)(i - 1)+ < s < 1, and < q < oo. Then 

v x :B^ q {dn)^B s p +Htt), 

Q-.B;-\dn)^B s p ^-\n), 

v x :B;jdii)^BZ*~\n), 



are well-defined, bounded operators. 
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2.2 Dirichlet data in L p (Q) 

We say a few words about the Dirichlet problem for the multi-dimensional Stokes 
system with L p boundary data which is not touched upon in [MMSlj . Following 
the paper by M. Dindos and V. Maz'ya |DMaj . we shall speak about both the Lame 
system (with the Poisson ratio a < 1/2) and the Stokes system (with a — 1/2). 

Let us consider a bounded domain fi in M. n , n > 3, and the system 

Am - Vtt = 0, div u + (1 - 2a)-n- = in n (2.6) 
complemented by the condition 

Tru = h e L p (dCl), (2.7) 

and the class of solutions is described by the inclusion 

u* e Lp(flfi). (2.8) 

The boundary values Tr u are understood in nontangential sense, that is in the sense 
of the limit 

Tru(x) = lim u(y), 

over a collection of interior nontangential cones r(:r) of some aperture and height and 
vertex at x € du) , and u* is the classical nontangential maximal function defined as 

u*(x) — sup for all x £ du. 

yer(x) 

One says that problem (|2.6|) . (|2.7[) is L p solvable if for all vector fields h £ L p (dft) 
there is a pair (u, tt) satisfying (|2.6[) - (|2.8p and moreover, for some C > independent 
of h, the estimate 

II u *IIl p (90) < C\\h\\ Lp{dn ) 

holds. Furthermore, the problem (|2.6[) . (|2.7|) is said to be solvable for continuous data 
if, for all h € C(9fi) the vector field u belongs to C(Tt) and the estimate 

l|w*llc(n) ^ C\\ h \\c(dn) 

holds. 

The L p solvability of the Dirichlet problem for the Stokes system is established by 
Z. Shen |Sh] for all p 6 (2 — e(n), oo] provided is a three-dimensional Lipschitz graph 
domain. In [DMa] . the problem (|2.6I) . (|2.7p is considered on domains in R™, n > 3, 
with isolated conical singularity (not necessary a Lipschitz graph) and the authors 
prove its solvability for all p € (2 — e(f2), oo] as well as its solvability in C(Q) for the 
data in C(dVl). This seems to be a strong indication that the range p 6 (2 — e(f2), oo] 
should hold also for Lipschitz graph domains. However, there is no such a result for 
n > 3. 

2.3 Lipschitz continuous solutions 

We cite a regularity result for solutions of the Dirichlet problem for the Stokes system 
in a plane convex domain obtained by V. Kozlov and V. Maz'ya in [KM3 : 

—Am + V7r = /, div u = in Cl 

(2.9) 

Tr m = on dn, 
where / € W^iVL) and (m,tt) e x L 2 (Q,). 
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Theorem 5 Let f2 be a bounded convex two-dimensional domain and let f € L q (Cl) 
for some q > 2. Then the velocity vector u € W^fil) admits the estimate 

l|Vu|| Loo( n)<C||/|k (n)) 

where C depends only on Q. 

A direct consequence of this result for the nonlinear Navier-Stokes system is as 
follows. 

Corollary 1 Let (u,ir) S W 2 1 (^) x L 2 (£l) solve the Dirichlet problem 

2 

—Au + V7r + UkdkU = /, div u = in f2, 

m = on dft, (2.10) 

where f £ W^ifl)- Let Q be a bounded convex two-dimensional domain and let 
f G L q (fl) for some q > 2. Then the velocity vector u £ W^i^l) belongs to the 
Lipschitz class C 0,1 (tt). 

A result of the same nature was obtained by V. Maz'ya in }Maz3| for solutions 
of the Neumann problem for the Poisson equation in arbitrary convex n-dimensional 
domain. Let £1 be a bounded convex domain in K™ and let W^(f2) stand for the 
Sobolev space of functions in L p (f2) with distributional derivatives of order I in W l v (il) . 
By L p ,_i_(f2) and W l j_(f2) one means the subpaces of functions in L p (f2) and Wj,(Q) 
subject to 

/ v dx = 0. 
Jn 

Theorem 6 Let f G L q with a certain q > n and let u be the unique function 
in W^lfl), also orthogonal to 1 in £2^) and satisfying the Neumann problem 

—Au = f in O, 

(2.11) 

^ = on on, 

au 

where v is the unit outward normal to dfl and the problem ^2.11\ is understood in the 
variational sense. Then the solution u £ j_(f2) of the problem (fJO jp satisfies the 
estimate 

l|V U || Loo( n) <C\\f\\ Lq[n) . (2.12) 

As a particular case of a result obtained by A. Cianchi and V. Maz'ya in [CM] 
for a class of nonlinear equations, one can replace the £ 9 -norm on the right-hand 
side by the Lorentz norm which is the best possible majorant among those 

formulated in terms of Lorentz spaces. Similar facts for Lame and Stokes systems 
with boundary conditions different from those of Dirichlet are unknown. 

3 Higher regularity of solutions 

3.1 Preliminaries 

In this section we give applications of Sobolev multipliers to the question of higher 
regularity in fractional Sobolev spaces of solutions to boundary value problems for 
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higher order elliptic equations in a Lipschitz domain. Since the sole Lipschitz graph 
property of £1 does not guarantee higher regularity of solutions, we are forced to select 
an appropriate subclass of Lipschitz domains whose description involves a space of 
multipliers. For domains of this subclass we develop a solvability and regularity 
theory analogous to the classical one for smooth domains. We also show that the 
chosen subclass of Lipschitz domains proves to be best possible in a certain sense. 

Let Q be a domain in R™ with compact closure. Throughout this section we assume 
that for any point of the boundary dtt there exists a neighbourhood U and a Lipschitz 
function ip such that 

U n fi = U n {z = (as, y) : x G R"" 1 , y > ip{x)}. (3.1) 

Let B l p ~ 1/p (W 1 - 1 ), I = 1,2,..., denote the completion of the space Cg°(H n - 1 ) in 
the norm 

( / / l^m-Mt) - iVpj-itiCaOHa! - tr +1 -^ Wi) 1/P + ||t»|| ip(B »-i). 
Vr«-i Je"- 1 7 

Replacing R" _1 by 90 one arrives at the definition of the space B p ~ 1 ' p (dCl). 

By the space of Sobolev multipliers M(W^(£1) —> W^(f2)) we mean the class of 
functions 7 such that ju G W^(0) for all u G The space M(W p ,l (0) -> W£(f2)) 

is endowed with the norm 

\\l\\M{w£(n)^wi p (n)) = sup{||7u|| w j(n) : < !}• (3-2) 

The notation MW^(O) is used for /i = I. Properties and applications of Sobolev 
multipliers are studied in |MS2] , 

We introduce the essential norm of a function 7 G M(Wp(Q) — > W p (Q.)): 

esa \h\\M(w£(a.)^wl(n)) = i*jf\h ~ T\\wg(a)^wi(n)i ( 3 - 3 ) 

where {T} is the set of all compact operators: Wp(Cl) — > Wp(O). 

Analytic two-sided and one-sided estimates for the norm (|3.2[) and the essential 
norm (|3.3|) can be found in |MS2j . 

3.2 Subclasses of Lipschitz graph domains 

Let H be a bounded Lipschitz graph domain. We introduce the class M p 1 ' p (I = 
2,3,...) of boundaries dfl, satisfying the following condition. For every point of d£l 
there exists an n-dimensional neighborhood in which <9f2 is specified (in a certain 
Cartesian coordinate system) by a function ip such that 

Vip e MB l p - 1 - 1/p (R n - 1 ). 

Furthermore, by definition, Mp is the class of bounded Lipschitz graph domains. 

We say that dfl belongs to the class Bp 1 ' p if dfl can be locally specified by a 
function p G B P ~ 1/p {W 1 ' 1 ). Since 

MB l - l - 1 lP{W l - x ) c B^ocT 1/P ( Rn " 1 )' 1 ^ 2 > 
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and C ' 1 ^™- 1 ) C B pA l' c p (W n ~ 1 ), it follows that any bounded domain D, with dtt G 



l ~ 1/p satisfies dn G B l p 1/p . 
According to Corollary 4.3.8 in IMS2], for p(l — 1) > n we have 



II^IIm^- 1 -^^"-!) ~ ™P_ 1 \\V<p\\ B l-i-Vp (Bl(x)) , 

where B r (x) is a ball in R n_1 with radius r and center x. Therefore, the classes 
Mp _1 ^ p and B P ~ X ' V coincide for p{l — 1) > n. 

For a bounded Lipschitz graph domain 17, by B l p ~ l ^ p (<917) we denote the space of 
traces on <917 of functions in W^(17). Taking into account an analogous fact for special 

Lipschitz domains of the class M p ~ x ^ p (see Subsection 9.4.3 in [M52]), we obtain that 
MB l p ~ 1 ^ p (dfl) is the space of traces of functions in MWp(fl). 

In our subsequent exposition the following additional condition on 17 plays an 
important role. 

We say that 517 belongs to the class M l J T 1 ^ p {8) if for every point O G 517 there 
exists a neighborhood U and a special Lipschitz domain G = {z — (x,y) : x G 
R n ~ l ,y > ip(x)} such that [/nfi = [/nG and 

Obviously, the boundaries in M P ~ 1 ^ P (S) belong to the class Mp _1//p and, therefore, 
to the class B l p ~ x ^ v '. 

The following assertion gives a local characterization of the class M p 1 ' p (6). In 
its statement we use the notion of the (p, j)-capacity of a compact set e in : 

C PJ (e) = infllMI^^ : u G Cg°(R n ^ 1 ), u > 1 on e}. (3.5) 

For various properties of this capacity see the books by V. Maz'ya |Maz2] and D.R. 
Adams and L.-I. Hedberg [AH] . 

Proposition 3 Letp(l — 1) < n. The class M p 1 ^ P (S) admits the following equivalent 
description. For any point O G <917 there exists a neighborhood U and a special 
Lipschitz domain G = {z — {x,y) : x G R n , y > <p(x)} such that U n 17 — U CiG and 

,. / \\Dt-x/ p (<p;B e )\\ L (e) \ j. 

( sup — + V^p Loo (B.) ) < c«y . (3.6) 

Here B £ is the ball with centre at O and radius e, c is a constant which depends only 
on I, p, n, and 



D^ 1/p ( v; B e )(x) = ( f \v 3 _M*)-v 3 -My)\ p \ x *l-2 +P ) 



dy \ 1 /p 
x-y\ 

Proposition [3] and properties of the capacity lead to the following sufficient condi- 
tion formulated in terms of the (n — l)-dimensional Lebesgue measure m n -\. 

Corollary 2 (i) If n > p{l — 1) and 

,. / IIA-i/p(y;S e )||i ( e) .. .. \ 

hm sup U ,, 's + \\\/ip r eg,) < cd , 
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then Oil £ M l p~ 1/p (5). 
(ii) If n = p(l — 1) and 



Now we present another test for the inclusion of a function into Mp 1 ^ p (6) involving 



lira (sup ||A-i/p(p;B B )|U p(e )|log(m n -ie)|&'- 1 )/P-|-||V V || ioo(1 ,. ) ) < cS , 

then dne M l ~ 1/p (5). 

Now we present 
the Besov space B™ . 

We say that the boundary of a Lipschitz graph domain f2 belongs to B l q ^ p (I — 
1,2,...,) if, for any point of dfl, there exists a neighborhood in which dft is specified 
in Cartesian coordinates by a function (p satisfying 

1 1 \ Vl - 1 ^ X + k ^- Vl - 1 ^\ qdX ) P/q \ h \2-n-p < °° ■ 

Corollary 3 Let p(l — 1) < n and let Q be a bounded Lipschitz graph domain with 
dVt G B l q ~p^ p , where 



l-i/p 

q e [ P (n - l)/(p(l - 1) - 1), oo] if p(l - 1) < n 



and 

g G (p, oo] if p(/ — 1) = n. 

Further, let dQ be locally defined in Cartesian coordinates by y — tp(x), where ip is a 
function with a Lipschitz constant less than cS. Then dfl G M p 1 ^ p (5). 

Setting q = oo in Corollary [3j one obtains the simple sufficient condition for the 
inclusion into Mp 1/,p (i5) formulated in terms of the modulus of continuity w/_i(t) of 

j^)'*<~ (3-7) 
Since B l ~% p C B p ~ 1/p , it follows that d3TTJ) is sufficient for dfl G B p ~ 1/p . 

3.3 General elliptic boundary value problem 

Consider either scalar or square matrix differential operators 

P(x,D x )= ]T a a (x)D%, Pj{x : D x )= £ a aj (x)D^ (3.8) 

|a|<2m |a|<fcj 

where x G f2, 1 < j < to, and D x — . . . , i _1 9 n ). We assume that the operators 

P, tr P\ , . . . , tr P m form an elliptic boundary value problem in every point O G d£l 
with respect to the hyperplane y = and that P is an elliptic operator in CI. 

The next result is proved essentially in the same manner as Theorem 14.3.1 in 
[MS2]. 

Theorem 7 (i) Suppose that for any neighbourhood U C M. n there exist operators 
P U (D Z ) = £ <%D„ Pf{D z ) = < D * 

|q|— 2m |a|— fcj 
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with constant coefficients such that {P u ; Tr PF} is the operator of an elliptic boundary 
value problem in the half-space {x € K™ _1 , y > 0}. 
(ii) Let 

\\a a (z) -OalUooCi/nn) < S, 

\a\—2m 

Yl eas ll°«lljif(wi-""(n)-»wJ-*»(n)) ^ S > ( 3 - 9 ) 

\a\<2m 

where I is integer, I > 2m, 1 < p < oo . The constant S here and elsewhere is supposed 
to be small in comparison with the norms of the inverse operators {P u ; Tr PF}^ 1 for 
all U. Further, let the coefficients a a j satisfy similar conditions with 2m replaced by 
kj. 

(Hi) Let the boundary d£l belong to the class M p 1 ^ p (8) if p(l — 1) < n or to the 
class Bp if p(l — 1) > n, in either case with 1 < p < oo. 
Then the operator 

m 

{P;TrPj} : W£(fi) -> W^ 2m (il) x J[ B^-^dQ) (3.10) 

i=i 

is Fredholm. In particular, for all u £ W p (il) the a priori estimate 

m 

Hlwj(n) < c (\\ Pu \\wl- 2m (n)+J2\\ TTP 3 u \\ B '-^-^p {dn) + Nk(n)) (3.11) 

holds; the last norm in the right-hand side can be omitted in the case of a unique 
solution. 

Note that even the roughest sufficient condition (|3.7p on the domain for (Hi) in 
Theorem [7] to hold is sharp. Let w be an increasing function in C[0, 1] such that 
w(0)=0, 

5 J u(t)-± + J u(t)j <cu(S), 

and 

It was shown in Section 4.4.3 of |MS2j that one can construct a function ip on R n 1 
such that 

(i) the continuity modulus of V;_i<p does not exceed cu with c = const; 

(ii) supp y> C Qi-n, where Qd = {x € R"^ 1 : < e?}; 

(*«) «p ^ Bp~ 1/p (R n ~ 1 ). 

Given 93, one can construct a bounded domain i7 in R n such that the Neumann 
problem 

Av — v = g in ft, dv/dv — h on Oft (3-12) 

with g e W l p - 2 (n) and h € B l p~ X ' 1/p (dQ.) may fail to be solvable in W l p (VL). 

Next we describe conditions on the coefficients of (13.8[) which are equivalent to 
those formulated in Theorem [7] and follow from the results in |MS2j . Ch. 7. In its 
formulation, we use the notion of capacity of a compact set in R" defined similarly to 
(ESJ). 
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Corollary 4 Conditions <\3.9^i in Theorem [7] can be stated as follows: 

(i) The coefficients a a with \a\ = 2m are in the class Wp~ 2m (£l) if p(l — 2m) > n 
and satisfy the inequality 

II ( \ U\\ , v ||Vj_2mOa|U„(e) , 

\\a a {z) ~ a a \\ Loo (unn) + ™ sup — , < <5 

e_>0 {eCfi:*am(e)< E } [^p,l-2m{^)\ /P 

if p(l — 2m) < n; 

(ii) The coefficients a a with \a\ < 2m are in the class Wp~ 2m (Q) if p(l — \a\) > n 
and satisfy the inequality 

-P Z^^Wp + SUP P 2m - H - n ' P \\ a MAB^)^<6 

e ^ QK {ecn:diam{e)<e} [^p,l- \a\ ( e )\ /P x£n,p<s ' 

if pQ ~ M) ^ n - Here B p (x) is a ball in M™ with radius p and center x. 

Various other sufficient conditions for (|3.9j) follow from the results in |MS2j . 

3.4 Dirichlet problem in terms of traces 

Let us first consider the Dirichlet problem as a particular case of the general boundary 
value problem dealt with in Subsection 3.3. We write the scalar or square matrix 
elliptic operator P in the form 

Pu= (-^) lal D a {A a p(z)D^u) . (3.13) 

\a\. |/3| <m 

Suppose that the coefficients A a p are in C l ~ m (n), l > m, and that the Garding 
inequality 

Re J A a pD a uTW^dz>c\\u\\ 2 wr(n) (3.14) 

^ \a\ — \/3\— m 

holds for u £ Co°(fi). Assume that f2 is a Lipschitz graph domain. 

We introduce a sufficiently small finite open covering {U} of f2 and a corresponding 
partition of unity {([/}• Let 

P jV = d j - 1 /dy j - x , j = l,...,m if U n dfl ^ 

and 

p jV = o if f/ n on = 0. 

The Dirichlet boundary conditions will be prescribed by the operators 

u 

We give a formulation of the Dirichlet problem. Let us look for a function u € 
WjJ(fi) such that 

Pu = / in ft, Tr Pju = fj on dQ, j = 1, . . . , m , (3.15) 

where / and /j are functions in Wp _2m (0) and B l p +1 ~ : '~ 1 ^ p (dfl) respectively. 

Here is the principal result relating the problem p. 151) borrowed from Subsection 
14.5.4 in [MS2]. 
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Theorem 8 Let any of the following conditions hold: 
(a) m = 1, p(l - 1) < n; 90 £ M^ 1/p (S); 
{fi) m = 1, P (l - 1) > n; 90 £ B l jT 1/p ; 

(7) m > 1, 90 € Mp 1 and 90 is locally defined by equations of the form 
y = <p(x), where tp is a function with a small Lipschitz constant (for p(l — 1) > n, this 
is equivalent to 90 £ B L p 1 ^ p ). 

Then the operator 

m 

{P;TtPj} : -> V^- 2m (0) x J[ B l p +1 ^- 1/p {dil) 

j'=i 

«s an isomorphism. 

An example in Subsection 14.6.1 of |MS2] shows that for p(7 — 1) < n and for 
m = 1 the condition 90 S M P 1//p (<5) in part (a) of Theorem [8] cannot be replaced by 
the assumption that 90 belongs to the class M l p 1 ^ p n C 1 ^ 1 . To be precise, a domain 

is constructed with 90 £ M^ 2 n C 1 for which the problem 

-Au = f in O, Tru = on 90 (3.16) 

is solvable in W^O) not for all / e £2(0). This means that the smallness of the 
seminorm || Vy|| a/b 1/2 (R"- 1 ) m ^ e definition of M% ($) is essential for the solvability 
of problem (13~T51) in Wf (0). 

The next assertion, which follows directly from the Implicit Function Theorem 
9.5.2 in [MS2] . shows that the condition 90 € B l v ljp with p(Z - 1) > n is necessary 
for the solvability of problem (|3 . 1 5[) in Wl(Q) for the operator P of higher than second 
order. 

Theorem 9 Let O be a bounded Lipschitz domain and let I be integer, I > 2m, 
p(l — 1) > n, 1 < p < 00, and m > 1. // t/iere exists a solution u £ W p (Vt) of the 
problem 

Pu = 0inO, Tru = 0, TiP 2 u=l, Tr Pjit = 0, j = 3, m , (3.17) 
i/ien 90 e P P ~ 1/p . 

Under the additional assumption 90 G C /_2,1 , the necessity of the inclusion 90 <E 
Bp~ 1/p for p(l - 1) < n is proved in Subsection 14.6.2 of [MS2] . 

Theorem 10 Let 90 oe m i/ie cZass C /_2,1 and let I be integer, I > 2m, p(l — 1) < n, 

1 < p < 00, and m > 1. // £/iere exists a solution u £ W p (Q) of problem §3. then 

90 £ B l - 1/P . 

A similar result for the second order operator P in Subsection 14.6.2 of IMS2j runs 
as follows. 

Theorem 11 Let I be integer, I > 2, 1 < p < 00, m — 1, and PI < 0. Let O be a 
domain with 90 £ C 1 and let the normal to 90 satisfy the Dini condition. If, for a 
nonpositive function f £ Cq°(0), there exists a solution u £ W^(0) of the problem 

Pu = f in 0, Trn = 0, (3.18) 

then 90 £ B l ~ 1/p . 



22 



Note that the convergence requirement (|3.7I) , the heaviest assumption on Q made 
in Corollary |31 is in a sense a sharp condition for solvability of the Dirichlet problems 
(|3.17l) and (|3.18[) in Wp(iY). The corresponding domain is constructed with the help 
of the same function ui as in the case of the Neumann problem (|3.12l) (see Example 
15.6.1 in [MS2] ), 

Example 15.5.2 from |MS2j shows that surfaces in the class M p 1 ^ p {5) with p(l — 
1) < n may have 

conic vertices if n > p(l — 1) 

and 

s-dimensional edges if s < n — p{l — 1). 

Suppose that for any point O <E d£l there exists a neighborhood U such that U (~1 f2 
is C°°-diffeomorphic to the domain 

R s x {{x,y) : y > tp{x s+1 ,.. . ,x n _i)}, < s < n- 2, 

i.e. the dimensions of boundary singularities are at most n — 1 — s. Then (|3.4[) is 
equivalent to 

and, in particular, it takes the form 

II^IIb'-^"^-!-*) ^ cS ' 

if n — s < p(l — 1) < n. In other words, dfl € M p ~ 1 ' p {8) if and only if (n — 1 — s)- 
dimensional domain {(a;, y) : y > tp(x s +i, . . . , a;„_i )} belongs to Mp" 1/p (c<5). 

3.5 Dirichlet problem in a variational formulation 

It turns out that for equations and systems of order higher than two, the formulation of 
the Dirichlet problem can be changed so that the solvability condition Q g M P 1 ' p (d) 
is replaced by the better one Q £ M P +1 ~ m ~ 1 ^ p (S). We comment on this referring to 
Section 14.5 of [MS2] . 

Let be open subset R n and let P be the operator (|3~T5)l . where A af} e C'~ m (fi), 
^ > to. Further, let the Garding inequality (|3.14p hold for u € Co°(0). 

We say that u € W^(O) is a variational solution of the Dirichlet problem if 

Pu = f, u-geWl(Q)nW^(fl) , (3.19) 

where / and g are given functions in the spaces Wp~ 2m (Q) and Wp(Q) respectively. 
By W~ k {Vi), k — 1,2,..., we mean the space of linear continuous functionals in 

We present an a priori estimate for solutions of problem (|3.19[) . 
Theorem 12 // 

(i) either p (I — to) < n, 1 < p < oo and dQ belongs to the class M P +1 m (6), 

or 

(ii) p(l — to) > n, 1 < p < oo, and dQ € B l p +1 m 1 ' p , 
then 

\\u\\wl{n) < c (\\ Pu \\w[,- 2m (n) + IMUi(n)) ( 3 - 20 ) 
for allue (w'nic)(i]). 
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Next we state two corollaries of (|3.20[) . 

Proposition 4 Let satisfy the conditions of Theorem ] 121 

(i) If the kernel of the operator 

p ■. (Wp n wp(n) -> w [ - 2m {vt) (3.21) 

is trivial, then the norm \\u\\L t (n) in !\3.20\ can be omitted. 

(ii) The kernel of the operator <\3.21\ is finite- dimensional and the range of this 
operator is closed. 

Proposition 5 Let f2 satisfy the conditions of Theorem ] 121 Further, let U and V be 
open bounded subsets ofR n , U C V and u G [W l p H W™)(tt). Then 

\\ u \\wi(Unn) < c (\\ Pu \\wl- 2m (vnfi) + IMUi(Vnn)) ■ 
Finally we give a theorem on the solvability of (I3.19[) . 

Let the Garding inequality (|3.14[) hold for all u € Cq°(CI). Then, as is well known, 
the equation Pu — f with / G M / 2 _m (ri) is uniquely solvable in W™(Sl). 

Theorem 13 Let dfl G M P +1 m for p(l — m) < n and let dfl belong to the class 
Bp for p(l — m) > n. 

00 if f e w l p - 2m {n) n w^ m (n), g e w*(ti) n w 2 TO (ft), 1 < P < oo, and if u e 

W^(Q) is such that Pu = f,u-g€ W^(fl), then u G W l p (0) and u - g G W™(Q). 
(ii) The problem \3.19§ has one and only one solution u G W p (Q). 

3.6 Strong solvability of the Dirichlet problem for the Stokes 
system 

Here we complement Section 2 by a solvability of the Dirichlet problem for the Stokes 
system in weighted Sobolev spaces of higher order. 

Let I be noninteger, I > 1. We use the notation M p (5) for the class of three- 
dimensional Lipschitz graph domains subject to 

II^IUs^- 1 ^) ^ 6 

for an arbitrary coordinate system on dtt, where S is a positive number. 

We conclude this section by stating strong solvability result for the Dirichlet prob- 
lem for the Stokes system 

Au — V7r = /, div u = g in Q 

(3.22) 

Tru = h on dQ. 

We assume that g and h satisfy the compatibility condition (|2.3p and use the space 
\ym,a introduced in Subsect. 1.3. The proof is essentially the same as that of Theorem 
15.1.2 in [MS2] . 

Theorem 14 Letp G (l,oo), a = 1 — {I} — 1/p, where I is noninteger, I > 1. Suppose 
that dtt G B l p for p(l - 1) > 2 and dQ G M l p (5) with some 5 = 5(p, I) for p(l - 1) < 2. 
Then, for every triple 

(f,g,h) g wf - x < a (n) x w^ a (n) x B l p (dn) 

there exists a unique solution (u,n) of the problem 22\ in Wp^ +1 ' a (Q) x Wp ,a (f2). 
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Note that under the conditions of the last theorem, the operator 

W^ +1 ' a {Q)xWl l] ' a (Cl)B(u,n) =^ 

(Am - Vvr, divu, Tiu) <E 1^-^(0) x W®' a (il) x B l p (dfl) 

is continuous. 

4 Asymptotic formulas for solutions to elliptic equa- 
tions 

4.1 Asymptotics of solutions near Lipschitz boundary 

Results of a new type were obtained by V. Kozlov and V. Maz'ya in |KM3j for solutions 
of the Dirichlet problem for higher order elliptic equations in Lipschitz graph domains. 
We mean an explicit description of the asymptotic behaviour of solutions near a point 
O of the Lipschitz boundary. 

Consider the special Lipschitz graph domain 

G = {x = (x\ x n ) G K" : x n > <p(x')}, 

where <p(0) — and ip has a small Lipschitz constant. The authors of |KM3) study 
solutions of arbitrary strongly elliptic equation of order 2m with constant complex- 
valued coefficients 

L(d x )u(x) = f(x) on B 3 n G (4.1) 

with zero Dirichlet data on (B 3 n dG)\0, where B p — {x : |x| < p} and d x = 
. . . ,d n ). It is assumed that the operator L has no lower-order terms and the 
coefficient in front of d^ m is equal to (—1)'". 

One of the results in KM3] is the existence of a solution U € W™(G) of the 
homogeneous equation (|4.1[) ) which admits the asymptotic representation 

U(x) = exp(- [ p(yi)d™E(y',0)dy' + O( [ x 2 (p)^)) 

K J\x\<\y'\<l J\x\ P ; 

(4.2) 

x f( Xn - V (x>)r + o(\xr+^( f K{ P )4r- e + 1))) • 

<*\x\ P 

Here e is a poisitive constant, 

x(p) = sup |V^(y')li 
\v'\<p 

and E is the Poisson solution of the equation L(d x )E(x) — in the upper half-space 
R™ which is positive homogeneous of degree m — n and is subject to the Dirichlet 
conditions on the hyperplane x n = 

d J n E = 0, 0<j<m-2, and d^E^Six'), 

where S is the Dirac function. 

It is also shown that a multiple of 11 is the main term in the asymptotic represen- 
tation of an arbitrary solution u if both u and / are subject to mild growth conditions 
near O. 
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Solutions with the weakest possible singularity at O are studied in [KM3] as well. 
The authors construct a solution IA of the homogeneous equation (|4.1I) which is subject 
to the asymptotic formula 



U{x) =exp( f V {y')d™E{y',Q)d y ' + 0{ f x 2 ( P )^-) 

X J \x\<\y'\<l J\x\ P 

x (E(x n - V (x')) + 0(\xr- n+1 -*( f x( P )4h + !)))■ 



(4.3) 



The asymptotic formulas (14.21) and (|4.3[) can be simplihed under additional condi- 
tions on x(p). Let, in particular, 

1 ..2,„\ d P 



x (p) — < oo. 

p 

Then, in the special case of the polyharmonic equation 

{-A) m u = on B 3 n G, 

any solution u satisfying |u(a;)| = 0(|x|" i_ ™ +1_e ) is subject to the following alterna- 
tives: either 

/ x (p(x')) m ( T(n/2) f . \ 

v I-™ 72 y|xi<i»'i<i \v \ ' 

or 

^)~C(^-^')rex P (-m^ / ^)|^)' 

v 71 ' J|a:|<| v '|<l \V I J 



4.2 Asymptotics of solutions to equations with discontinuous 
coefficients near a smooth boundary 

Proofs of the just mentioned results in |KM3j rely upon the papers [KM1 , KM2 
on the asymptotic formulas for solutions to the Dirichlet problem for arbitrary even 
order 2m strongly elliptic equations of divergence form near a point O at the smooth 
boundary. It is required only that the coefficients of the principal part of the operator 
have small oscillation near this point, and the coefficients in lower order terms are 
allowed to have singularities at the boundary. 

We say a few words on the proof of asymptotic formulas in [KMl] , IKM2] . The el- 
liptic equation is transformed to a first-order evolution system with the matrix whose 
entries are partial differential operators on the hemisphere with time dependent coeffi- 
cients. Thus, the question of asymptotics of solutions to the original Dirichlet problem 
is reduced to the study of the long-time behaviour of solutions of the evolution system 
just mentioned. The structure of the operator matrix in the system is rather compli- 
cated, because it has been obtained from a higher order partial differential equation in 
the variational form. Moreover, the study of this system is aggreviated by the scant- 
iness of information about the behaviour of the operator matrics at infinity. This 
difficulty is overcome by a right choice of function spaces, characterizing the solutions 
and the right-hand side of the evolution system by certain seminorms depending on 
time. To obtain an asymptotic formula for the solution, the authors apply a par- 
ticular spectral splitting of the system into one-dimensional and infinite-dimensional 
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parts. A general asymptotic theory of differential equations with operator coefficients 
in Banach spaces which is the basis of [KM1 , KM2 is developed in [KM] . 

As an illustration, we describe a corollary of the main result in |KMlj concerning 
second order equations. Consider the uniformly elliptic equation 

- div (A(x) Vtt(x)) = f(x) in fl (4.4) 

complemented by the Dirichlet condition 

u = on dQ, (4.5) 

where il is a domain in R n with smooth boundary. We assume that elements of 
the n x n-matrix A{x) are measurable and bounded complex-valued functions. One 
consideres a solution u with a finite Dirichlet integral and suppose, for simplicity, that 
/ = in a certain (^-neighborhood £1$ — {% € f2 : |x| < 6} of the origin. Further, it is 
assumed that there exists a constant symmetric matrix A with positive definite real 
part such that the function 

cr(r) := sup \\A(x) - A\\ 

is sufficiently small for r < 6. 
We introduce the function 

((A(x) - A)u, v) - niA-^Aix) - A)v,x){u,x){A~ 1 x,x)- 1 



n(x) 



\S n ~ 1 \(detA) 1 / 2 {A- 1 x,x) n /' 2 



where \S n 1 1 is the Lebesgue measure of the unit sphere in R™ , (z,Q = ^iCi + - • -+ z nCn 
and v is the interior unit normal at a point O on the boundary of tt. (For the notation 
(detA) 1 / 2 and (A~ 1 x,x) n / 2 see 0, Sect. 6.2). 

The following asymptotic formula for an arbitrary solution of (|4.4|) . (|4.5j) with 
finite energy integral is a special case of the main theorem in [KM1 : 
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u(x) = Cexp(- I n(y)dy + 0{[ <r( P ) 2 —)) 

V Jflg\fl\B\ J\<H\ P ' 

(dist^, an) + o{\ x \ 2 -- J s <t(p)-^)) + o(\ x \ 2 -n, (4.6) 



x 



where C — const and e is a small positive number. 

Using (|4.6p , it is an easy matter to derive sharp two-sided estimate for the Holder 
exponent of u at the origin. Another direct application of (14.61) is the following 
criterion. Under the condition 



5 a (pf^P <O0 (4.7) 
o P 



all solutions u are Lipschitz at the origin if and only if 



hO 



lim inf / ReKdx>-oo. (4.8) 



n s \n T 



Needless to say, this new one-sided restriction (|4.8I) is weaker than the classical Dini 
condition at the origin. The complementary assumption (I4.7|) appeared previously in 
several papers dealing with other problems of the boundary behaviour of solutions to 
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equation (|4.4|) (see the articles by E. Fabes, D. Jerison, and C. Kenig |F JKj . by B.E. 
Dahlberg |Dah3j . by C. Kenig [Ken] et at). 

Let v be a solution of the equation (|4.4j) complemented by the Dirichlet condition 

v = on dn\{0}, O € dil, 

which has an infinite energy integral and the least possible singularity. We state 
a particular case of the main theorem in |KMlj which is the following asymptotic 
representation for v: 

v(x)=Cexp([ K{y)dy + 0(f a(pf^)) 
vn,\(i|»i J \A P ' 

where C — const and e is a small positive number. 

In general, theorems proved in [KMlj and |KM2) provide asymptotic formulas 
similar to (|4.6j) and (|4.9[) for solutions of the Dirichlet problem for the higher order 
uniformly elliptic equation with complex- valued coefficients 

(-d x ) a (£ a0 (x)dPu(x)) =f(x) on B+, 

0<|a|,|/3|<m 

where B+ = R'| n B s , = {x = (x' , x n ) e R n : x n > 0} and B s = {x e K n : \x\ < 
5}. The only a priori assumption on the coefficients C a p is smallness of the function 

\ a \ = \f3\ =m \a+/3\<2m 

where x € Bj and L a p are constants. 

4.3 Corollaries of the asymptotic formulas in Section 4.1 

The last section in [KM3 concerns, in particular, solutions to the Dirichlet problem 
for elliptic equations of order 2m with constant coefficients in plane domains with a 
small Lipschitz constant of the boundary as well as arbitrary bounded plane convex 
domans fl. Let fi be a bounded domain in K 2 . Consider a strongly elliptic operator 
with constant coefficients 

L(d x )= ]T a k d k x dl m -\ 

0<fc<2m 

and denote by w a weak solution to the Dirichlet problem 

L{d x )w = /, we w?(n). (4.io) 

If / € W2 _m (r2), this problem is uniquely solvable, let us assume that 
CinB 2 6 = {{xi,x 2 ) : x 2 > <p(xi), \x\ < 2S }, 

where if is a Lipschitz function on [—2Sq, —25q] and <p(0) = 0. Note that one does 
not require the convexity of tp. 

The next result concerning solutions to problem (14. 10)) , which stems from (|4.2I) , is 
as follows. 
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Theorem 15 Suppose that the Lipschitz norm of ip on [— 28o, — 2t5o] does not exceed 
a certain constant depending on the coefficients of L. Let f be equal to zero in fJni^cS- 
Then, for all S € (0, So), x G f2 D B$ and k = 1, . . . , m — 1, 

\V k w(x)\<cA(2S)\x\ m - k (4.11) 

xexp (_ a f m^t£l dp + h / 5 ma x|^(i)| 2 ^Y 
v J\x\ P J\x\ |t|<p P ' 

Here 

A(S) = 8- 1 -™ |MU a (nn*,) (4.12) 

and the notation 

a = ^ E (tf-C). 

l<fc<m 

is used, where £^ and are roots o/ the polinomial L(l,£) u>z£/i 

positive and negative imaginary parts, respectively. This value of a is best possible. 
By b and c positive constants depending only on m and the coefficients of L are 
denoted. 

Note that for the operator A m one has = ±i, which implies a — —m/ir. 

The next assertion is a consequence of Theorem 1151 when the function ip is convex. 

Theorem 16 Suppose that the function ip describing the domain Q near the point O 
is non-negative and convex, and |<//(±2<5)| does not exceed a sufficiently small constant 
lo depending on m and the coefficients of L(d x ). Furthermore, let f be zero in QHB26 
and let w be a solution of l\4-.l(ty , which is extended by zero outside fl. Then 

\\V m w\\ Loe(Bs) <cA(85), (4.13) 

where S < So/ 8 and A{5) is given by \4- Ity - 

One of the main results obtained in [KM3J concerns the Green function Gl of the 
Dirichlet problem for the operator L with real coefficients. 

Theorem 17 Let Q, be an arbitrary bounded convex domain in R 2 and let the coeffi- 
cients of L be real. Then, for all x, y in £1, 

£ \d«d?G L (x,y)\<C\x-y\- 2 , (4.14) 

\oi\= |/3|— m 

where C is a positive constant depending on f2. 

The case of complex coefficients is more complicated. 

Theorem 18 Let L be an arbitrary strongly elliptic operator with complex coeffi- 
cients. Suppose that Q is a bounded convex domain such that the jumps of all angles 
between the exterior normal vector to dfl and the x-axis be smaller than a constant 
depending on m and the coefficients of L(d x ). Then, for all x, y in Q, estimate 
holds. 

Theorem [T"7l implies the following regularity result. 
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Corollary 5 Let Q be an arbitrary bounded convex domain in M 2 and let the coeffi- 
cients of L be real. Then the solution w of problem \4-.10\ with f £ Wq m (Cl), q > 2, 
satisfies 

E \\D a w\\ LUQ) <C\\f\\ (4.15) 

a|<m 

where the constant C depends on £1, m, q, and the coefficients of L(d x ). 

Generally, this assertion does not hold for operators with complex coefficients. 
More precisely if there exists an angle vertex on the boundary of a convex domain 
Q, one can construct a second order strongly elliptic operator L(d x ) with complex 
coefficients such that the Dirichlet problem (|4.10p with / e C(fl) has a solution with 
unbounded gradient (see [KMR] . Sect. 8.4.3). By Theorem H8l the statement of 
Corollary [5] for L with complex coefficients holds if the jumps of the normal vector 
are either absent or small. 



4.4 Classical asymptotics of solutions near a point of the do- 
main 

Now, we present some results borrowed from Sect. 14.9 of the book by V. Kozlov 
and V. Maz'ya [KMj which are devoted to the asymptotic behaviour of solutions to 
elliptic equations near an interior point O of the domain. Here a modified Dini-type 
condition on the coefficients is introduced which guarantees the preservation of the 
asymptotics of solutions to the main part of the equation with coefficients frozen at 
O. 

Let 

P(D X )= E P« D % 

\a\—2m 

and let G denote the Green matrix of this operator, i.e. the solution of the system 

P{D X ) G(x) = It S(x) inM" (4.16) 

where /; is the I x I identity matrix and S is the Dirac function. It is well-known (see 
F. John's book [Jo]) that G admits the representation 



G(x) 



r 2m u Q(uj) if 2m > n, n odd, or 2m < n 

R(x) logr + r 2m ~ n S(oj) if 2m > n, n even, 



where Q and S are smooth matrix-functions on the unit sphere in M. n and R is a 
homogeneous polynomial matrix of degree 2m — n. 

Let us consider the elliptic operator 

Q(x,D x )= £ q a {x)D a x 

\a\ <2m 

with measurable coefficients in the punctured ball i? ro \{0}. We introduce the function 
S(r)=sup{ £ \q a (x)~ Pa \+ £ \x\ 2m -^\q a (x)\], 

r \a\—2m \a\<2m 

where K r = {x E W l : e~ 1 r < \x\ < r}, and assume that S(r) does not exceed a small 
positive constant. We shall also use the notation 

\a\<2m 
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We formulate three theorems on the asymptotic representation as x — > for solu- 

F 

Q(x, D x )u = on B ro \{0} (4.17) 



tions u e W 2 2 r oc (^o\{0}) of 



satisfying 

NUfW = 0(r fe + 5 ) (4.18) 

with some 5 > and integer A:. 



Theorem 19 Lei 2m < n and 



f r ° dr 
/ S(r)|bgrr- x -< 
Jo r 



(4.19) 



where J is a positive integer. 

(i) lfk>0, then 

u{x) = c a x a + v(x), (4.20) 

|a|=fe+l 

where c a — const and 

||«|| W | m( x r .)=o(r fe+1 |logr| 1 -^). (4.21) 

(ii) If k < 2m — n — 1, i/ien 

u(ss)= 51 C a D«G(x) +v(x), (4.22) 

|a|— 2m — n — fc— 1 

where C a = const and G is the Green matrix introduced by \4- 

(iii) If k = 2m — n, then 

u(x) = const + v(x), (4.23) 

where 

\\v\\ wim(Kr) = odlogrl 1 ^). 

The asymptotics (I4.23P can be made more precise under the assumption that the 
operator Q contains no derivatives of order \a\ < s, that is 

Q(x,D x )= ]T q a {x)D a x . 

s<\a\<2rn 

The formula (|4.23p can be replaced by 

u(x) = 5^ c a x a + v(x), 

\a\<s 

where 

ll«llw|«(Jf P ) =o(r s |logr| 1 - 7 ). 
Theorem 20 Let n be odd, 2m > n and let S be subject to <\4- Ity - Then 

u(x)= J2 c <x x ™ + J2 C/3D^G(x)+v(x), (4.24) 

|a|=fc+l |/9|=2m-n-fe-l 

where c a and Cp are constants and v satisfies <\4-21]) . If either k < — 1 or k > 2m — n, 
then the first or the second sum in Q^.^P should be omitted. 
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Theorem 21 Let n be even, 2m > n. 

(i) If k < —2 and <\4- 19\ holds, then u satisfies <\4- 14\ with v subject to \^.21\ . 

(ii) If k > 2m — n and <\4- 19\ holds, then u satisfies \4-2Q[ with v subject to l\4-2ty . 
(Hi) Let — 1 < k < 2m — n — 1 and let \4.19\ be valid with 7 > 2. Then u is 

represented by \4-%4\ with 

\\u\\ wlm{Kr) =o(r k+l \\ogrr-<). 

4.5 Asymptotics of solutions of the second order equation with 
square-Dini coefficients 

The asymptotic behaviour of solutions near the isolated point O of the domain was 
recently considered by V. Maz'ya and R. McOwen [MMlj . |MM2j for the case of the 
second order elliptic operator in nondevergence form 

C(x,D x )u(x)— aij(x)didjU. (4-25) 

l<i,j<n 

It is assumed that the coefficients have modulus of continuity ui satisfying the square- 
Dini condition 

/ Lu 2 (t)— < 00. (4.26) 
Jo t 

If the coefficients are real, then, without loss of generality one can put -(0) = &y. 
An important role is played by the function 

' {r) '\^\L i J t '^ ) - n{J W L )w (4 - 27) 

where A 2 stands for the matrix aij(z), (,) is the inner product in R", and s is a 
sufficiently small positive number. 

If the coefficients are subject to the usual Dini condition 

/ ui(t)— < 00, (4.28) 
Jo t 

then, obviously, there exists a finite limit of I(r) as r — > 0, but (|4.28|) is not necessary 
for the existence of this limit. In general, under (|4.26|) . I(r) may be unbounded as 
r — s- 0, but, clearly, for every A > there exists C\ such that 

|J(r)| < A|logr| + C X for < r < e. (4.29) 
The results in jMMlj and |MM2j are formulated in terms of the L B -means 

M p (w,r) := (-/ \w{x)\ p dx) 1/P . (4.30) 

^J r<\x\<2r ' 

The main theorem in [MMlj can be stated as follows. 
Theorem 22 Suppose that 

\a-ij(x) — Sij\ < u(\x\) as x O, 
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where uj satisfies \^.2b\ . For p e (l,oo) and e > sufficiently small, there exists a 
weak solution Z € L p ^ oc (B e ) of 

d i d j (a ij (x)Z)=0 (4.31) 

l<z,j<n 

satisfying 

Z(x) = e- I ^(l + C(x)) 1 (4.32) 
where the remainder term £ satisfies 



M,,i^.r) -. niwxuirUW! witli a(r! := / —f^ dt - (4.33) 



o 

Moreover, if u € £ P ,; O e(i?£\{0}) is a weafc solution of 

^ d l d j (a ij (x) u) = 

l<z,j<n 

m i? e subject to the growth condition 

M p (u,r) < cr 2 - n+eo , 
where Sq > 0, then there exists a constant C , depending on u, such that 

u(x) = C Z(x) +w(x), (4.34) 
where the remainder term w satisfies 

M p (w,r) < cr 1 ' 51 

for < r < e and any E\ > 0. 

Combining (|4.34l) and (I4.32[) . one obtains the asymptotic representation 

u(x) = ce- I ^\x\ n (l + o(l)) asx^O. (4.35) 

Now let us turn to the behaviour of a solution to the homogeneous nondivergence 
equation with an isolated singularity at x = O. We assume that n > 2. The results 
obtained in [MM2 imply the existence of the solution 

Z{x) ~ ^—r e /(|x|) as x -> O. (4.36) 
n — 2 

The behaviour of /(|a;|) as x — > O not only controls the leading asymptotics of Z{x) 
but whether we can solve the equation 

jC(x,D x )F(x) = S(x). 

There are three important cases to consider. 

1. 1(0) = lim iY|x|) exists and is finite. 

In this case Z[x) may be scaled by a constant factor to make it asymptotic to the 
fundamental solution for the Laplacian. In fact, the distributional equation 

-C(x,D x )Z(x) = CS(x) (4.37) 
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can be solved to find 



C= IS^le 1 ^. 



(4.38) 



2. I(\x\) — > -oo as x -» O. 

We see that 

Z{x) = o(|x| 2 ~") as x O, 

and we can solve (|4.37|) to find C = 0. Thus, in this case we obtain the interesting 
corollary that 

Cu = in B £ (4.39) 
admits a solution u = Z that is quite singular at x = O: 

\Z(x)\ > C x \x\ 2 - n+x 

for every A > 0. In particular, local regularity of solutions of the homogeneous 
equation (|4.39[) does not hold. 

3. I{\x\) — > oo as x — > O. 

Now we find that 

Z(x)\x\ n ~ 2 ->• oo as i ^ O, 

so this solution grows more rapidly than the fundamental solution for the Laplacian. 
Although Z still satisfies (|4.37p . we can no longer find C. 
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